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ABSTRACT 

Calibration is one of the long-standing problems in optical interferometric mea- 
surements, particularly with long baselines which demand stars with angular sizes on 
the milliarc second scale and no detectable companions. While systems of calibra- 
tors have been generally established for the near-infrared in the bright source regime 
(K< 3 mag), modern large interferometers are sensitive to significantly fainter magni- 
tudes. We aim at providing a list of sources found unresolved from direct observations 
with high angular resolution and dynamic range, which can be used to choose interfer- 
ometric calibrators. To this purpose, we have used a large number of lunar occultations 
recorded with the ISAAC instrument at the VLT to select sources found to be unre- 
solved and without close companions. An algorithm has been used to determine the 
limiting angular resolution achieved for each source, taking into account a noise model 
built from occulted and unocculted portions of the light curves. We have obtained up- 
per limits on the angular sizes of 556 sources, with magnitudes ranging from K s ft4 
to 10, with a median of 7.2 mag. The upper limits on possible undetected companions 
(within ft 0''5) range from K s ft8 to 13, with a median of 1 1.5 mag. One-third of the 
sources have angular sizes < 1, and two-thirds < 2 milliarcseconds. This list of unre- 
solved sources matches well the capabilities of current large interferometric facilities. 
We also provide available cross-identifications, magnitudes, spectral types, and other 
auxiliary information. A fraction of the sources are found to be potentially variable. 
The list covers parts of the Galactic Bulge and in particular the vicinity of the Galactic 
Center, where extinction is very significant and traditional lists of calibrators are often 
insufficient. 

Subject headings: Techniques: high angular resolution - Occultations -Stars: binaries: 
general - Stars: fundamental parameters 
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1. Introduction 



Lunar occultations (LO) can efficiently yield high angular resolution with short observations 
using relatively simple instrumentation. As such, they have been widely employed to measure 
stellar angular diameters and detect small-sep aration binary sourc es , providing the bulk of such 



information for several decades. The catalog (|Richichi et al. 



2005) listed several hundreds of LO 



results. However, the introduction of large format detectors in the visual and near-IR has eroded 
the opportunities at most telescopes for fast time resolution as required by LO. At the same time 
long baseline interferometetry (LBI) has progressively filled the requests for the same kind of 
fundamental observations, without the LO limitations on the choice of the sources. 

LO have recently been re-introduced as a productive method to obtain high angular resolution 
measurements that may not be collected efficiently in any other way. By using fast integrations 
on subwindows of array detectors in the near-IR, LO have been observed for several hundreds 
of so urces at the ESQ Ve ry Large Telescope (VLT) using a minimal amount of telescope time 



(e.g. 



Richichi et al. 



20121 . and references therein). LO targets are selected randomly according to 



the lunar apparent orbit, and naturally the majority of them turn out to be unresolved. Nevertheless, 
these unresolved sources are also of interest, in particular for what concerns interferometry. 
Indeed, the angular resolution achieved in these recent LO observations at the VLT is of order 1 
milliarc second (mas) and thus comparable to the typical LBI resolution on a large interferometer, 
while the sensitivity is typically several magnitudes better (Km 12 mag). 

Therefore, sources which are found to be unresolved by LO are well suited to be used as 
calibrators by near-IR LBI. Considering that in our experience about 10% of the randomly- selected 
LO sour ces turn out to b e resolved or binary, and that so-called "bad" LBI calibrators are a 



concern (Hummel 



2008), such a selection of clean calibrators is indeed desirable. In this work, 



'Based on observations made with ESO telescopes at Paranal Observatory 
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we list 556 sources found to be unresolved from our LO observations at the VLT, and thus highly 
homogeneous in terms of data and analysis. For each, we provide an upper limit on the angular 
size as well as a sensitivity limit on the possible presence of undetected companions. The sources 
are mainly located in the Galactic Bulge, with several in the vicinity of the Galactic Center where 
visual extinction can be extreme and therefore it might be more difficult to identify suitable LBI 
calibrators. 



2. Observations and data analysis 



The o bservations which form the basis of this work were obtained using the ISAAC 



instrument (Moorwood et al. 



1998) at the ESO VLT. The instrument was installed at the 8.2-m 



UT1 Antu telescope until May 2009, and thereafter at the twin UT3 Melipal telescope. The 
majority of the observations were carried out in service mode following a filler strategy to utilize 
small amounts o f time available bet ween other programs or due to unfavorable atmospheric 



conditions (e.g., 



Richichi et al 



in visitor mode (e.g., 



20101). A significant amount of observations was also carried out 



Richichi etalJl201lh 



Each observation consisted of several thousands of frames in a small sub-window with short 
sampling (also equal to integration) times. Typical values were 7000 frames of 32x32-pixel 
(4"7 x 4 "7) at 3.2 ms, although different combinations exist across our data set. A broad-band K s 
filter was generally employed, except in the case of very bright sources for which a narrow-band 
filter centered at 2.07 jum was used. The events were generally disappearances, although a number 
of reappearances were also observed. As expected from the nature of the data set, the observing 
conditions covered a very wide range, with median values of seeing, airmass, lunar phase being 
0"88, 1.30 and 63%, respectively. The final quality of the light curves is only marginally affected 
by atmospheric conditions, however other important factors are the image quality (LO at the VLT 
are observed with the primary mirror active optics off) and the intensity of the lunar background. 
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These can vary significantly from one LO event to another and may be responsible for a large 
scatter in signal-to-noise ratio (SNR) among light curves of sources with similar magnitude. 

The data cubes were converted to light c urves using a mask e xtraction, and these in turn 



were analyzed u sing bot 



m ethod (CAL, 



in 



i a m odel-dependent (IRichichi et al. 



Richichi 



Richichi et al 



199q) and a model-independent 



1989). Details on the instrumentation and the method can be found 



(1201 1|) and references therein. Large samples of LO data inevitably include 



occasional cases of light curves with slightly distorted diffraction patterns, most likely as a result 
of irregularities in the lunar limb. This appears to be statistically more frequent in observations 
with a large telescope w hich samples a wider section of the limb and which also provides 



an increased sensitivity. 



Richichi et al. 



(|2012|) discussed the nature and the effect of such 
irregularities, and suggested correction algorithms that we employed for some light curves in our 
sample. 

An important aspect of our data reduction was the systematic assessment of the limiting 
angular resolution achieved for each of the light curves. This was based on an algorithm which 
evaluates the quality of the fit over a range of angular diameters, taking into account also factors 
such as the finite number of dat a points and their intrinsic numerical accuracy. Details are 



provided in 



Richichi et al. 



(I1996|) . This quantity, which we denote with U , can be considered 
equivalent to an upper limit on the angular size of the source. In Fig. Q] we show an example of the 
data and best fit by a point-source for the case of a high SNR light curve, while Fig. [2]provides a 
similar exmaple for a low SNR case. 

For each light curve, we combine the SNR and the K magnitude to derive the detection 
limit K lim for a hypothetical companion. In a LO event the noise level might not be symmetric 
around the time of occultation. Also, depending on the circumstances and geometry, detections 
of companions below the SNR=1 level can be possible. Therefore, the Ki; m quantity is not a 
strict limit but rather an indication. We note that only restricted portions of the light curves, 
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2000 2200 

Relative Time [ms] 

Fig. 1 . — The top panel shows the occultation light curve (dots) and the point-source fit (solid 
line) for 2MASS 17383557-2223159, a K=6.0mag source without known cross-identifications. 
The SNR of the fit is 148 and the limiting angular resolution is <p u = 0.50 mas. The lower panel 
shows the fit residuals, enlarged for clarity. 
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corresponding to angular extensions of « 0"5, were analyzed. For our scopes, the definition of 
unresolved source is only meaningful within this restricted range. Companions further away than 
this limit have been noticed in some cases or might have been present without being detected, 
however we do not concern ourselves with them. 



3. Source selection 

A total of 915 successful LO events were logged in our database from March 2006 until June 
2012. A significant fraction of these observations yielded positive results in areas of binary stars, 
stellar diameters, and sources with extended circumstellar emission are described in other papers 
already published or in preparation. Here, we concern ourselves only with the sources effectively 
found to be unresolved. Several tens of sources already foun d to be resolved or binary during our 



20121 and references 



observations at the VLT have already been published (e.g. iRichichi et al 
therein). 

The first step in their selection was the rejection of all light curves showing binary or multiple 
stars, within the angular range mentioned in Sect. [2] This included borderline cases in which 
the CAL algorithm revealed a secondary peak in the brightness profile, generally close to the 
noise level, but for which no significant improvement in the quality of the fit could be reached 
with a double-star model. We then rejected all stars with a resolved angular diameter or an 
otherwise extended apperance, including those for which the <p u algorithm yielded an indication 
of being resolved even though no diameter could be fitted. We also rejected light curves with 
significant distortions of the fringe pattern due to presumed limb irregularities, e.g., exceeding by 
several times the noise level over several fri nges. However , mod erate distortions were generally 



satisfactorily treated with the algorithms of 



Richichi et al. 



(|2012|) and the corresponding light 



curves were retained. Finally, a minimum SNR of 10 was imposed. 
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As a result, a total of 556 sources were selected and are listed in Tabled] available online. 
We used the 2MASS Catalogue for our predictions, and further identifications are extracted from 
the Simbad database. We also provide our own ID number that combines the ESO Period number 
with the entry in our database. This ID also marks with an asterisk sources observed with a narrow 
band filter. Further we list coordinates, BVJHK photometry, spectrum, number of references and 
source type, as available from Simbad. The next fields in the Table are the SNR, K lim and <p u as 
described in Sect. 2 and computed by us. Finally, we flag sources known or likely to be variable. 
By the flag S we denote sources which belong to categories such as Miras, emission-line stars, 
RR Lyr or other kind of variables. By the flag F we denote stars for which our measured counts 
differed from the counts expected on the basis of the 2MASS K s value by more than 0.6 mag (see 
also Fig.©. 



4. Discussion 

The distribution in magnitude for the sources in our list is shown in Fig. [3j We used the 
nominal K s values from 2MASS, although our measurements have shown at times significantly 
different values as discussed later. They range from K s =3.7 to 10.2, with a median of 7.2 mag. 
The figure also shows the distribution of the limiting magnitudes Kn m as introduced in Sect|2] The 
values range from K s =8.2 to 12.8, with a median of 1 1.5 mag. 

The distribution of the upper limits on the angular sizes <f> u is shown in Fig. 0] The values 
range from 0.25 to 8.1 mas, with a median of 1.4 mas. A total of 192 sources, or about one third 
of the sample, have a limiting angular resolution of 1.0 mas. The (fi u < 2.0 mas limit is satisfied 
for 371 sources. Fig. [5] shows another view of the cp u distribution, this time as a function of 
the K s magnitude. For this, we converted the measured counts into an observed K s magnitude, 
using the ISAAC exposure time calculator (see also Fig. [8]). The shaded area in Fig. |5] covers 
the range between the minimum and the median <p u values over one-magnitude bins. The scatter 
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in </> u is a reflection of different SNR encountered for a same source brightness depending on 
observing conditions, as mentioned in Sect.[2l However, when just the best 50% of the sources are 
considered in each magnitude bin, the LO performance appears to be quite consistent and shows 
that the 1 mas limit is attained up to K?a 7 mag. 

It is interesting to compare our upper limits U with the predicted angular diameters. Several 
emp irical relation s have been derived, to estimate a star's diameter from its colors. A recent one 



by 



Kervella et al. 



(2004) is based on V and K photometry and can be applied to dwarf stars in 
the spectral range AO to M2, and sub-dwarf stars between AO and K0. Unfortunately, we face 
the complication that several of our stars suffer from considerable reddening, have little spectral 
information, and in some cases appear to be variable - as explained in the following. Nevertheless, 
we have endeavoured to determine the estimated diameter for a subset of 90 stars within the above 
constraints of photometry and spectral class, without taking into consideration possible interstellar 
reddening. The result is that in all cases the U is larger than, or at the most approximately similar 
to, the estimated diameter. Only in one case we found some discrepancy (our upper limit being 
about 25% larger than the estimated diameter). This is IRAS 14434-2055, a source with over 4 
magnitudes of V-K color, which is likely not intrinsically photosperic. We conclude that, as far as 
it is possible to check, our upper limits are well consistent with the expected stellar diameters. 

The majority of the sources have no known counterpart, with just under 50% having any 
identification other than the 2MASS name. About 30% have at least one bibliographical entry, 
and only 13% have at least two. Spectral types are known for only 15% of the sources, and V 
magnitudes for only one out of three. This is a reflection of many of the sources being in the 
Galactic Bulge, see Fig. [6l where interstellar extinction is significant or even extreme. As a result, 
the sources in our list are also mainly red or very red in color, as shown in Fig.[7J While generally 
consistent with a giant sequence, the near-IR colors show that many of the sources are affected by 
10 or even 20 magnitudes of visual extinction. 
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In fact, we note that 10 sources are within 1° of the Galactic Center, and 45 are within 5°. 
This highlights a potential use of our list, namely to provide reference stars for LBI near-IR 
observations in a very interesting area of the sky, for whi ch calibrators selec ted by other methods 



such as the modelling of the spectral energy distribution (|Cohen et al. 



1999) are very hard to find. 



Fig. [8] shows the average unocculted counts within our extraction masks for all sources, as 
a function of the nominal 2MASS K s magnitudes. It can be observed that in general the counts 
follow very closely the theoretical ISAAC performance. However, in a significant number of 
cases deviations are observed both in the sense of the observations showing brighter or fainter 
magnitudes than 2MASS. Small variations can be attributed to airmass (which we did not account 
for), atmospheric conditions and to the mask extraction process. However, larger deviations are 
likely indicative of intrinsic variability. In Fig. [8] we have used 0.6 mag as a threshold, leading to 
about 9% of the sources appearing to have varied their brightness since the 2MASS measurement. 



The problem of finding suitable calibrator stars in LBI has been long-standing. 



3) 



Cohen et al 



1999) produced an extensive compilation of predicted angular diameters, based on modelling 



of the spectrophotometric properties of bright non- variable stars. This has often been used as 
a source of calibrators, and vari ations more suitab le for the specific application of near-IR LBI 



have also been produced (e.g., 



Merand et al. 



2006J). More recently, as extensive databases of 



observations have become available from facilities such as the Palomar Testbed Interferometer and 
the ESO Very Large Telescope Interferometer, attempts have been made t o derive self-consisten t 



systems of ca 



Richichi et al 



ibrato rs with measured, as opposed to predicted, diameters (|van Belle et al 



2009) 



2008 



While the choice of suitable calibrators is largely eased by the above mentioned works, a few 
limitations are still encountered: for example, the density of calibrators is less than satisfactory 
in some areas of the sky, and they are usually limi ted to rather br i ght s tars, typically K< 3 mag. 



Tools for selecting fainter calibrators are available 



Bonneau et al. 



(|201ll) : however they are again 
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based on indirect methods and they may require interpolations for missing photometry. The 
quoted photometric accuracy of 0.1 mag is generally acceptable to estimate angular diameters, and 
large relative errors are not a problem in the faint magnitude (small size) regime. However, such 
accuracy does not allow us to ascertain in advance, for example, whether a specific source has a 
faint companion. On the other hand, experience from our LO observations has shown that, with a 
resolution of « 1 mas and a dynamic range of AKss 5 mag, about 10% of the field stars appear to 
have companions. 

Interferometric observations are generally quite time-consuming both in acquisition and 
analysis of the data, and it is often only after the observing run has been concluded that problems 
connected with a calibrator, such as a previously unknown companion, can be recognized. In 
this sense, our list of unresolved sources can prove useful. Admittedly, we provide only upper 
limits rather than angular diameters, however for baselines which are not too long this can be 
sufficient. For example, on a 50 m baseline in K band, the difference in squared visibility between 
a 1 mas and a point-like source is 3%, comparable to the typical accuracy of an interferometric 
measurement. Of course, LBI observations repeated on different baselines can detect deviations in 
the point-like nature of a source with better accuracy. The case of an unknown faint companion, 
however, might be harder to recognize. 



5. Conclusions 

We have used a database of about 900 near-IR lunar occultation light curves obtained with 
the ISAAC instrument at the VLT, to identify a set of sources that can be classified as unresolved. 
This set comprises 556 stars spanning a range of magnitudes from Ks=4 to 10, which have been 
analyzed in a highly homogeneous context in terms of instrument performance and data analysis. 
We are able to set upper limits on their angular sizes, typically of 1 to 2 milliarcseconds, and 
on the presence of nearby companions, with an average dynamic range of 4.5 mag, reaching a 
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typical limiting magnitude Ki^ll to 12. The sources are distributed along the Zodiacal belt, and 
especially clustered in the Galactic Bulge. They exhibit mostly red or very red colors due to 
interstellar extinction. A significant number are in close proximity to the Galactic Center. As a 
result, only a small fraction of the sources in our list are well studied, with the majority lacking 
extensive classification and bibliographical references. 

We identify a key application of this list as a database of calibrators for long-baseline 
interferometry, especially in the faint source regime of modern large facilities and for targets in 
the general direction of the Galactic Center where other sources of reliable calibrators are scarce. 

The observations which form the basis of this work would not have been possible without 
the support and dedication of the ESO staff, both in Garching and in Chile, whom we thank 
sincerely. OF acknowledges financial support from MICINN through a Juan de la Cierva 
fellowship and from MCYT-SEPCYT Plan Nacional I+D+I AYA#2008-01225. MM acknowledges 
co-funding under the Marie Curie Actions of the European Commission (FP7-COFUND) through 
the ESO Fellowship program. This research made use of the Simbad database, operated at the 
CDS, Strasbourg, France, and of data products from the Two Micron All Sky Survey (2MASS), 
which is a joint project of the University of Massachusetts and the Infrared Processing and 
Analysis Center/California Institute of Technology, funded by the National Aeronautics and Space 
Administration and the National Science Foundation. 
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2000 2200 2400 

Relative Time [ms] 

Fig. 2.— Same as Fig.Q] for 2MASS 06251004+2302141, a K=8.2mag source coincident with 
TYC 1879-97-1. The SNR of the fit is 16 and the limiting angular resolution is <fi u = 3.05 mas. 



-14- 



"i 1 1 1 1 1 1 r 



350 
300 
250 
200 
150 
100 
50 


3 4 5 6 7 8 9 10 11 12 13 14 

2MASS K s (mag) 

Fig. 3. — Distribution of 2MASS K s magnitudes (red, filled) and limiting magnitudes K lim (green, 
striped) for the sources in our list. 
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Fig. 4. — Distribution of limiting angular resolution, <f) u , for the sources in our list. 
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Fig. 5. — Distribution of limiting angular resolution </> u for the sources in our list, as a function of 
the observed magnitude. The shaded area is delimited at the bottom and at the top by the minimum 
and the median, respectively, of the U values over one-magnitude bins. 
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Fig. 6 



— Sky distribution of the sources in our list. The Galactic Center is marked with a cross. 
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Fig. 7. — 2MASS color-color diagram for the sources in our list. The 
unreddened giant (red) and dwarf (blue) stars, respectively, acc ording to 



The arrow is the extinction vector for A v =10mag, according to 



ines are the loci of the 



Bessell & Brett! (1988). 



Rieke et al.l ( 



1985) 
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2MASS K s mag 

Fig. 8. — Expected (lines) and measured (dots) flux values for the 2MASS K magnitudes of the 
sources in our list. The expected values are based on the ISAAC exposure time calculator (solid 
line for broad-band K s , and dashed line for narrow-band 2.07 /im filters). The circled dots denote 
sources which we flag as variable, see text. 
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Table 1 . List of sources found unresolved by the lunar occultation technique 



Object 


Our ID 


Simbad ID 


q(J2000) 


S(J2000) B V 


J 


H 


K Sp #Ref 


Type 


SNR 


Klim 


4>u 


Var 


17383141-2812445 


P76/gc097 


DENIS-PJ173831.3-281244 


17:38:31.4 


-28:12:44.5 


8.88 


7.47 


6.78 1 


IR 


70.8 


11.4 


1.50 




17385533-2836105 


P76/gcl37 


2MASSJ17385533-2836105 


17:38:55.3 


-28:36:10.5 


9.96 


8.33 


7.34 1 


IR 


44.8 


11.5 


1.80 




17385407-2833227 


P76/gcl47 


DENIS-PJ173854.0-283322 


17:38:54.1 


-28:33:22.7 


9.80 


8.03 


7.13 1 


IR 


47.5 


11.3 


1.20 




17394022-2820124 


P76/gcl99 


2MASSJ 1 7394022-28201 24 


17:39:40.2 


-28:20:12.4 


9.81 


8.12 


7.12 1 


IR 


84.3 


11.9 


3.70 


1 


17394602-2822089 


P76/gc207 


2MASSJ 1 7394602-2822089 


17:39:46.0 


-28:22:08.9 


9.72 


7.87 


7.05 1 


IR 


92.8 


12.0 


2.35 




17405112-2817130 


P76/gc260 


2MASSJ 1 7405 112-2817130 


17:40:51.1 


-28:17:13.0 19.10 


7.07 


5.51 


4.87 1 


IR 


181.7 


10.5 


0.60 


F 1 


17411415-2818051 


P76/gc284 


2MASSJ1741 1415-2818051 


17:41:14.2 


-28:18:05.1 18.50 


8.88 


7.49 


6.87 1 


IR 


100.4 


11.9 


0.75 


F 


17413435-2829225 


P76/gc352 


ISOGAL-PJ174134.6-282922 


17:41:34.3 


-28:29:22.5 


6.96 


5.52 


4.85 5 


Star,IR 


241.7 


10.8 


0.50 


F 


17425497-2836593 


P76/gc426 


2MASSJ17425497-2836593 


17:42:55.0 


-28:36:59.3 


10.49 


7.90 


6.72 1 


IR 


92.8 


11.6 


0.80 




17430901-2837039 


P76/gc437 


2MASSJ 1 743090 1-2837039 


17:43:09.0 


-28:37:03.9 


10.59 


8.16 


6.99 1 


IR 


85.7 


11.8 


0.75 





Note. — 



Table[T]is published in its entirety in the electronic edition of the Astrophysical Journal Supplement Series. A portion is shown here for guidance regarding its form and content. 



